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Abstract. We examine by extensive computer simulations the self-diffusion of 
anisotropic star like particles in crowded two-dimensional solutions. We investigate the 
implications of the area coverage fraction cj) of the crowders and the crowder-crowder 
adhesion properties on the regime of transient anomalous diffusion. We systematically 
compute the mean squared displacement (MSD) of the particles, their time averaged 
MSD, as well as the effective diffusion coefficient. The diffusion appears ergodic in 
the limit of long traces, such that the time averaged MSD converges towards the 
ensemble averaged MSD and features a small residual amplitude spread of the time 
averaged MSD from individual trajectories. At intermediate time scales we quantify 
the anomalous diffusion in the system. Also, we show that the translational—but not 
rotational—diffusivity of the particles D is a non-monotonic function of the attraction 
strength between them. Both diffusion coefficients decrease as D{(j)) ^ (1 — 
with the area fraction (p occupied by the crowders. Our results might be applicable to 
rationalising the experimental observations of non-Brownian diffusion for a number of 
standard macromolecular crowders used in vitro to mimic the cytoplasmic conditions 
of living cells. 


1. Introduction 

Over the recent years, deviations from the standard Brownian diffusion law [1] have 
been observed in a broad range of systems [aEllllElEllTllH]. Depending on the physics 
of the system under consideration, various theoretical models are used to describe these 
deviations [21 El H El El m E]. Such anomalous diffusion is typically characterised by 
the power-law growth of the mean squared displacement (MSD) of particles with time 

{r\t)) ~ (1) 

We distinguish subdiffusion for 0 < /? < 1 and superdiffusion for 1 < (3. Subdiffusion 
is an abundant phenomenon for passive motion in the world of live biological cells 
[H El El El E] • In the biological context subdiffusion was observed for particles ranging 
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from small proteins [S] via messenger RNA molecules HD] in the cell cytoplasm to large 
chromosomal loci and telomeres in the nucleus m to sub-micron virus particles [12] as 
well as lipid granules PI- The features of anomalous diffusion depend on the energy 
landscape and the physico-chemical interactions in the system of particles m US]. 
The advances of modern single particle tracking experiments Uni [ig [IT] provide a 
wealth of high resolution experimental data to quantitatively compare the microscopic 
mechanisms of non-Brownian diffusion with known theoretical models. The latter 
include, inter aha, the continuous time random walk [IHlIiniEolEIlESI or the equivalent 
formulation in terms of fractional diffusion equations [3l[23], fractional Brownian motion 
[21], heterogeneous diffusion processes PI. scaled Brownian motion [26l ET] 128], as well 
as the fractional Langevin equation related to the viscoelasticity of the environment 
[291130]. 

The cytoplasm of biological cells is a superdense PI fluid consisting of proteins, 
nucleic acids, membranous structures, cellular machinery components, semiflexible 
hlaments, etc. [SHEa Eg El]. This macromolecular crowding (MMC) reaches volume 
occupancies oi (j) > 30% [35]. In addition, the cytoskeletal meshwork [36] of eukaryotic 
cells impedes the diffusion of larger entities in cells, in particular, near the cell’s plasma 
membrane. The cytoplasm in addition is highly heterogeneous both in prokaryotic and 
eukaryotic cells [371 EH |39]. The anomalous diffusion of cell-related phenomena may 
represent a blend of more than one theoretical model representing the quality of the 
diffusion on different length and timescales [H El El 001 011 021 03l 011 05] . 

A number of experimental [371 06], theoretical 071 , and simulation [Ml 0E1 091 
EniEIlE2lE3lElE5lEniE7] studies in recent years were devoted to tackling various 
aspects of particle diffusion in crowded environments. From the simulation perspective, 
for instance, the studies of tracer diffusion in non-inert [58] , heterogeneously distributed 
and poly-disperse [59], restrictively mobile [M], squishy m and anisotropic [60l [61] 
obstacles were performed. Despite the progress of analytical theories of crowded 
solutions some important diffusive characteristics can only be studied quantitatively 
by computer simulations. This is particularly true for crowders of the non-trivial of the 
Mercedes-Benz® star like particles considered in the current paper (Fig. [1]). 

We here use computer simulations to unravel the implications of the particle shape 
and squishiness as well as the crowding fraction on the translational {D) and rotational 
{Dr) particle diffusivities in highly crowded solutions. Our main target is to gain 
insight into the physical behaviour of non-spherical crowders relevant for the situation 
in vitro where soft non-spherical and often non-inert crowders such as globular PEG 
and branched dextran polymers are routinely used to mimic the effects of MMC in 
living cells. Another important experimental example is the diffusivity of anisotropic 
lysozyme-like proteins studied by Brownian Dynamics simulations in crowded media 
[62]. It was demonstrated that—particularly in heavily crowded solutions—not only a 
transient subdiffusion of the protein centre of mass exists, but diffusion becomes also 
progressively anisotropic. This anisotropy of the translational diffusion pronounced on 
short-to-intermediate times disappears in the long time limit. The long time diffusivity 
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Figure 1. (A) Mercedes-Benz® star shaped crowder, with the centre monomer in red 
and flexible arms in blue. Typical conformations of crowders for (B) purely repulsive 
and (C) attractive interactions of strength ca = at crowder fraction (f> = 0.12. 

Video files illustrating the dynamics of the stars at €a = 0,1, 2kBT are provided in the 
Supplementary Material. 


values were shown to drop drastically with the protein concentration [62]. Moreover, 
the reduction of Dr for Y-shaped proteins such as IgG y-Globulin (molecular weight 
MW~ 155 kDa) was shown to be stronger than for more spherical proteins such as 
Bovine serum albumin (MW~66 kDa). These experimental observations based on 
fluorescence correlation spectroscopy measurements are supported by all atom Brownian 
Dynamics simulations [62]. The inclusion of hydrodynamic interactions revealed an 
additional reduction of Dr of proteins [63]. The reader is also referred to the simulation 
study of Ref. [M] in which the self diffusion of star like polymers in the presence of 
hydrodynamic interactions [65] was examined in detail. 

The paper is organised as follows. In section [2] we introduce our simulation model, 
the physical observable we are interested in, and some details on the data analysis 
algorithms. We present the main hndings of simulations in section |3l In section ID the 
implications of our results for some cellular systems are discussed. 

2. Simulation model and observables 

We implement our computer code developed to simulate the particle diffusion of crowded 
solutions in which all particles are explicitly treated [52] (53] |5l]. Here, we consider a 
two-dimensional system of Mercedes-Benz® star shaped crowders, each consisting of 
four discs of diameter a connected by elastic springs, see Fig. [T]A.. The elastic potential 
between the midpoint of the molecule and the centres of the outer monomers is 

Ucir) = ^ks{r-r^y, ( 2 ) 

where Vc is the equilibrium distance and kg the spring constant. We also connect the 
outer monomers with springs of the force constant kg, namely, 

Uo{r) = ^kg{r - Vof, 


(3) 
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Figure 2. Translational and rotational time averaged MSD of star like crowders for 
varying strength of the inter-particle attraction strength ca- For the time averaged 
MSD S'^{A) only the x components 5^{A) are shown—the y components show identical 
features. The insets show the translational and rotational particle diffusivities in the 
long time limit. Parameters: cj) = 0.15, T = 2 x 10®, the average (i5^(A)) is computed 
over iV = 40 traces. 


to mimic the softness of our triangular star like crowders. The equilibrium distances 
and constants are set to = l.Sa, Vg = 1.5\/3(j, and kg = IOO/cbT/ct^. The interaction 
between all beads is described by the 6-12 Lennard-Jones potential 


f^Lj(’^,’^cut) = 4e 


-(7 + 

r / 


12 


©(’^cut -r) -hC'(rcut)- 


(4) 


Here Q{x) is Heaviside step function and C(rcut) is a constant setting > rcut) = 0. 

For a purely repulsive potential the standard cutoff distance rcut = 2^/®cr is used with 
the potential strength e = ksT. For attractive interactions we set rcut = 2(J with varying 
adhesion strength e = between the monomers. This attraction acts between all the 
monomers of the stars. We use periodic boundary conditions within a square box of 
area L^. The packing fraction of N crowders in the system is dehned as 0 = NA/L"^, 
where A = 47 r((j/ 2 )^ is the total area of the four monomers and N ~ 10^ is a typical 
number of stars used in our simulations. In most scenarios below the system size is 
L = 40a and the total simulated trace length is ~ 4 x 10® of elementary time steps. 

The dynamics of the two-dimensional position rj(f) of the ith monomer disc 
interacting with the other monomer discs is described by the Langevin equation 

~ + Uoffij) + f/Lj(rp)] + (5) 

j 

Here ^{t) represents Gaussian white noise with zero mean {$,{t)) = 0 and correlator 
{^(t) ■ ^{t')) = 4'-ykBT6{t — t'), where ks is the Boltzmann constant, 7 is the friction 
coefficient, and T the absolute temperature. In the following, we use a and ksT 
as the basic units of length and energy, respectively. We simulate the system with 





























Self-subdijjusion in solutions of star-shaped crowders 


5 


the Verlet velocity algorithm with elementary time step At = 0.005 for the total 
time T. The physical time scale in these simulations is the standard combination 
6t = (J^Jm/{kBT) 0.36 ns [66], if we set the monomer diameter to a = 6 nm 
and its mass to the average mass of cytoplasmatic crowders, MW~ 68 kDa [STl 167] . 

We track the positions of the centre monomers of all the crowder stars and their 
orientation with respect to the x-axis, denoted as 9i. From the trajectory of the ith 
crowder we calculate the time averaged translational {6'^) and rotational MSDs 

as [7| 

+ A) -Xiit)]"^ + [yi{t + A) -yi{t)f}dt 
= + ^ly ( 6 ) 

and 

y.(A,r) = +A) - ()At)rdt. (7) 

Here A is the lag time along the trace. In addition to the individual MSDs 6‘^{A) we 
compute the corresponding averages over the set of N individual trajectories, 

1 ^ 

(i^(A)) = -j;i?(A), (8) 

' i=l 

as well as their amplitude spread around this mean value. 

The diffusion is called ergodic if the ensemble and time averaged MSDs coincide 
in the limit A/T —)■ 0 and if the spread of 5^ around the mean approaches the delta 
function in this limit [7]. A more accurate description of ergodicity can be achieved 
based on the so-called ergodicity breaking parameter EB. The latter is dehned as the 
variance of the distribution of the dimensionless variable f = 5^/ whose precise 

behaviour as a function of the lag time and the various model parameters, however, is 
beyond the computational scope of the current investigation. 

3. Results 

In Fig. [2] we present the behaviour of the translational and rotational MSDs for varying 
inter-particle attraction strength eA at crowder packing fraction 0 = 0.15. The initial 
crowder diffusion is ballistic, stemming from the simulation of inertial particles, see also 
Ref. [58]. At intermediate time scales of A ~ 0.1... 10 we observe a non-monotonic 
behaviour of the time averaged MSD that we ascribe to the events of the first collision 
of a given crowder molecule with another crowder. We quantify the variation of the 
local scaling exponent in Fig. [6f3 below. In the long lag time limit the translational 
and rotational MSDs grow linearly with A reflecting the Brownian behaviour of the 
crowder particles. In this limit the diffusion is ergodic, as we demonstrate in Fig. |3l 
This statement is not necessarily trivial: in many weakly non-ergodic systems the time 
averaged MSD turns out to be a linear function of the lag time A while the ensemble 
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Figure 3. Individual time averaged MSD traces and their dependence on the 
trajectory length T, plotted for the parameters of Fig.[2]and ea = 2kBT. The ensemble 
averaged MSD is the bold black line. The time averaged MSD trajectories become more 
reproducible for longer trace length T. 


averaged MSD scales as a power-law or logarithmically in time t. This phenomenon 
was observed in various experiments [T^ 021 03l 01] and explained in terms of various 
stochastic processes [H |25l |26l |68l |69l [TOl [711 ESI [73]. Fig. 0] demonstrates both that 
to very good approximation ergodicity in the sense of the equality {x^{A)) = 6‘^{A) 
and that a very small amplitude scatter around the mean (^^(^A)) exists and thus the 
time averages are reproducible quantities. We furthermore detail the dependence of the 
particle diffusivity in this Brownian limit versus the attraction strength and the Filing 
fraction in Figs. 0]and[71 respectively. 

Let us be more specihc. Fig. [3] illustrates the time averaged MSD for different 
lengths of the time series of the star diffusion as well as the superimposed ensemble 
averaged MSD shown as the bold black line. As can be seen from the figure, the 
amplitude scatter of single traces 5^ around their mean remains small along the entire 
trajectory except when A ~ T, as expected. This growing spread as A ~ T is a 
standard feature of even canonical Brownian motion appearing due to progressively 
poorer statistics when taking the time average [7]. More importantly we observe that 
the amplitude spread of time averaged MSD at a fixed lag time A decreases as the 
length T of the time traces increases. This property is ubiquitous for ergodic diffusion 
processes [7]. We note that the magnitude of the amplitude scatter that we observe 
for (5^ for moderately adhering Mercedes-Benz® stars are similar to that of a tracer 
in a network of sticky spherical obstacles, compare Fig. [3] above and Fig. 7 as well as 
the explanations in the text of Ref. [58]. Computing the magnitude of the mean time 
averaged MSD for varying trace length T we observe that its magnitude stays nearly 
unchanged with T (Fig. [3]). In the short lag time regime the ballistic scaling is visible. 
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Figure 4. Average Brownian diffusivity of crowders measured along the cc-direction 
[D = Dx here) versus their mutual attraction strength eA/iksT), plotted for the 
parameters of Fig. and varying crowding fractions f. The error bars are included 
and are often smaller than the symbol size. 
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Figure 5. Effective potential (HU) between the crowders for different attraction 
strengths eA/iksT), plotted for the parameters of Fig. [2] using the same colouring 
scheme. 

Given these observations our process is ergodic and thus fundamentally different from 
other anomalously diffusive systems such as those described by continuous time random 
walks or heterogeneous diffusion processes [25]. For the latter a pronounced scatter of 
the time averaged MSD trajectories around their mean and a clear dependence of the 
amplitude of ((5^(A) ) on T at fixed A exist, that is, the system ages [U [7ll 175] . 


The particle diffusivities are defined as 




2A 

( 9 ) 



<1 

CM 

( 10 ) 


obtained in the long time limit A 1. Fig. 0] shows the values of D and Dr obtained 
from a linear ht of the translational and rotational time averaged MSDs in the range 
A = 10^ ... 10^. We find that while the rotational diffusivity Dr decreases monotonically, 
the translational diffusivity remarkably exhibits a shallow yet significant maximum at 
~ IksT. This systematic trend persists for the variation of the crowder fraction 
in a quite broad range (Fig. 0]). This implies that the self-diffusion of our star like 
crowders can be facilitated by a weak inter-particle attraction. This is one of the main 
conclusions of this study. One can rationalise this trend in the self-diffusion in terms of 
the concept of the effective crowder size that decreases for moderate attraction strengths 
eA ~ 1/c_bT. Fig. 0] illustrates that for progressively stronger star-star attraction their 
mutual diffusivity decreases eventually to zero due to aggregate formation, see also Fig. [5] 
and its discussion below. 
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We also detect a progressive aggregation of crowders at relatively large values of 
the crowder-crowder attraction strength e^, as demonstrated in Fig.[Tp and in the video 
hies in the Supplementary Material. This is a well known phenomenon, for instance, in 
the glass transitions of dense suspensions of sticky hard spheres [76]. Accordingly, the 
decrease of the average diffusivity D of crowders as plotted in Fig. [Demerges due to the 
averaging over an ensemble of particles that perform individual random motions. This 
average takes into account both particles forming transient aggregates as well as free 
particles. Roughly speaking the average diffusivity drops inversely proportionally to 
the number of particles in the cluster. The fraction of particles clustering in these 
aggregates increases with the mutual attraction strength. The average diffusivity 
therefore progressively decreases with due to a larger fraction of particles in the 
transient aggregates. At large attraction strength the majority of particles belong to 
big clusters. 

At a hxed cohesiveness of our Mercedes-Benz® stars vicinal crowders create 
a rough energy landscape for the self-diffusion and the hopping of a given crowder 
particle. As the MMC fraction increases, the binding events give rise to more a 
prolonged particle aggregation and reduced self-diffusivity. Above a critical MMC 
fraction the barrier height exceeds the thermal energy unit thus increasing the 
lifetime of crowder aggregates signihcantly. For stronger star-star attraction the 
formation of essentially permanent aggregates sets in for less crowded systems leading 
to an inhomogeneous, phase-separated spatial distribution, see Fig. [TjC. 

A similar non-monotonicity of the translational diffusivity D at similar strengths 
of the particle-crowder attraction was found in Ref. [5T| for the tracer diffusion in dense 
suspensions of spherical Brownian particles. While we here detect that the attraction 
strength yielding the highest value of D is a function of the crowding fraction (f of 
the stars for the spherical particles, the stickiness facilitating the particle diffusivity 
was almost 0-independent in Ref. ISH. The non-monotonic D^sa) dependence was 
interpreted in Ref. [51] in terms of the roughness of the free energy landscape for the 
tracer diffusion using the concept of the chemical potential. Interestingly, the tracer 
diffusivity was also non-monotonic in 0 in a static regular array of sticky obstacles, as 
quantihed in Ref. [58]. 

We checked the universality of the observed dependencies for D^ca) and T*(0) also 
for spherical particles. Namely, we simulated just a single monomer of our Mercedes- 
Benz® stars with the given adhesive properties. The diffusivity was indeed found to 
reveal a maximum at ~ 0.5... 1 /c_bT (not shown), indicating some universality of 
this a priori counter-intuitive faster diffusion for a weak inter-particle attraction |51j . 
Note also that for a polymer chain diffusing in an array of sticky obstacles a weak 
chain-obstacle attraction can also substantially enhance the polymer diffusivity [5I1177]. 

To rationalise the observed behaviour of D{eA) we calculate in Fig. [5] the potential 
of the mean force between two crowders as 


F(r) = -ksT \og[p{r)]. 


( 11 ) 
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In this reconstructed free energy p{r) is the average radial distribution function of the 
centre monomer of the crowders in the steady state long time limit. As the mutual 
attraction strength increases we observe that the potential well at the separation 
r ^ a becomes deeper, see the hrst well in Fig. [5l Concurrently, the distance at 
which F{r) sharply increases becomes shorter for larger e^. For a stronger star-star 
attraction the crowders feature a more organised appearance, resulting in measurable 
oscillations of p{r) and F{r), as evidenced in Fig. [5l These trends indicate that the 
effective crowder radius gets smaller with increasing e^, and at an optimal value the 
crowders approach one another more closely yet without sticking. This in turn might 
result in a faster average diffusivity H at ~ e2’^ as we indeed observe. An effective 
reduction of the crowder size at optimal attraction strength is one important cause— 
albeit possibly not the only one—for this facilitated diffusion. In the current system, 
the equilibrium distance of the outer monomers from the cenral monomer is reduced 
by about 2% for the inter-monomer attraction strength of 2 /cbT. Even higher values 
of eA give rise to the formation of large clusters of crowders, see the Supplementary 
Material. As shown in Fig. IHthe corresponding diffusivity of an average particle drops 
dramatically. 

In Fig. [6l\ we show the translational and rotational MSD for varying packing 
fraction of crowders 0. As expected—from a linear ht to the long time time averaged 
MSD—the diffusivity is a monotonically decreasing function of 0, as evidenced by Fig. [71 
For more crowded systems the tracer diffusion gets more obstructed and the magnitude 
of the corresponding mean time averaged MSD decreases. To elucidate these effects 
further we evaluate from the time averaged MSD traces of Fig. [6l4 for the translational 
motion the local diffusion exponent [7117H] 

dlog[(¥{A))) 
dlog{t) 

For the rotational motion the exponent I3r{t) is defined analogously. We observe a 
ballistic regime with 0 = 2 in the particle diffusion at short times, (Fig. [6fl). This 
ballistic regime is followed by a decrease and further increase of the scaling exponent at 
A ~ 1. These non-monotonic trends are also clearly visible from the behaviour of the 
time averaged MSD traces themselves as a function of the lag time A, see Fig. [6]A. We 
hnd that the variations of the scaling exponent for translational and rotational motions 
of the star like crowders appear correlated, indicating a coupling of these diffusion modes 
[T6] . In the plots for the scaling exponent 0(t) in Fig. [6|B the signihcant spike-like signal 
at A ~ 1 is interpreted as an effect of the first collision of particles and the resulting 
onset of an effective confinement. We note that even in effective one-particle theories 
pronounced oscillations occur at the crossover point between initial ballistic and the 
overdamped regime [79l [80] . 

With increasing crowder fraction 0 we also observe a more pronounced range of 
anomalous diffusion for lag times of the order of A ~ 1... 100. This range appears 
strongly correlated between rotational and translational particle motion, as shown in 
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(A) 


(B) 



Figure 6. A: Translational time averaged MSD of the central star monomer and 
rotational time averaged MSD of the star polymer. B: local scaling exponent /3(t) of 
the time averaged MSDs computed for varying packing fractions f. In B, in the limit 
of short times a linear sampling of data points was chosen for the left panel and a 
logarithmic sampling for the right panel. Parameters: eA = IksT, T = 2 x 10®, and 
N = 10. 


Fig. [6]B. For rotational diffusion the scaling exponent drops practically to zero for times 
longer than those of the initial ballistic growth, and the corresponding mean time 
averaged MSD trace exhibits a short plateau (Fig. EB)- In the long lag time limit 

the exponent becomes Brownian ~ 1. Such a transient subdiffusion was observed for 
a number of systems see also the Introduction. Especially in dense colloidal 

systems close to the glass transition 0 = 0* this subdiffusion is accompanied by an 
exponential growth of the solution viscosity t] = ri{(j)) which is divergent at 0 —)■ 0* |81j . 
The colloidal glasses also exhibit progressive particle localisation effects as discussed in 
Refs. [HU [82]. 

Remarkably, the relative variation of the translational and rotational diffusivities 
with the crowding fraction of stars is quite similar. For comparison, we plot in Fig. [7| 
the theoretical prediction for dense suspensions of hard spheres 


D{0) 


1 - 


(13) 


with the critical packing fraction for our system of 4>*{eA = ^ksT) ^ 0.52. Above 
this value 0* both translational and rotational diffusivities of the crowders essentially 
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Figure 7. Normalised translational (measured along the x-direction) and rotational 
diffusivity as function of the crowding fraction (j) computed for ca = ^ksT, together 
with the asymptote of Eq. (1131) with D{(j) = 0) = Dq. 


vanish. At this critical crowding fraction the inter-particle attraction becomes so 
strong that the self-diffusion is almost completely localised and the motion of particles 
corresponds more to a very restricted wiggling and jiggling. As expected, when the 
star-star interactions become stronger aggregate formation sets in for less crowded 
systems, and thus the critical value <p* is diminished (not shown). Albeit this theory 
in Refs. |H3l [HI] is developed for three dimensional suspensions in the presence of 
hydrodynamic interactions, it agrees remarkably well with our results, as shown in 
Fig. [71 The reader is also referred to Ref. |85] for experimental data of the crowding 
dependent diffusivity of colloidal particles and alternative theoretical predictions for the 
diffusivity We note that Ref. [86| suggest exponential rather than power-law 

forms for the particle diffusivity in crowded solutions. 

4. Discussion and conclusions 

We performed extensive computer simulations and theoretical data analysis of the 
diffusion of crowders with a branched structure. A simple example of such spiky 
but responsive crowders in two dimensions are deformable Mercedes-Benz® like stars 
employed here. Their outer monomers are inter-connected by an elastic potential 
bestowing upon it a certain degree of responsiveness—an important characteristics for 
many polymeric crowders m- We also incorporated in the simulations an inter-particle 
attraction strength which represents another realistic feature of solutions of non-ideal 
crowders in vitro. 

We found that the diffusion of our Mercedes-Benz® star like crowders is ergodic 
and, within accuracy. Brownian in the long time limit. We examined the behaviour 
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of the ensemble averaged MSD and the time averaged MSDs of the crowders in a 
wide range of the crowder fraction cj) and the inter-crowder attraction strength e^. 
As a function of the crowding fraction we demonstrated that both translational D 
and rotational Dj. diffusivities follow the analytical decrease (IT^ of -D(0) predicted 
for suspensions of hard spheres. The dependence of the star-star attraction strength is 
more remarkable. Namely, the translational diffusivity shows a weak yet systematic non¬ 
monotonic dependence on for the solutions at all crowding fractions studied herein. 
The rotational diffusivity, in contrast, is a monotonically decreasing function of the 
inter-particle attraction strength e^. Thus, a relatively weak inter-monomer attraction 
can facilitate the lateral diffusion and also induce a certain degree of clustering and 
spatial heterogeneities in crowded solutions of non-inert particles. These effects will 
impact the diffusion of a tracer particle in crowded solutions—such as those of PEG, 
dextran, or Ficoll—used in vitro to mimic the crowded conditions in living cells [171119]. 
In addition to the proof of the ergodic long lag time diffusion shown in Fig. [3] and the 
transient subdiffusion regime of our star like crowders in Fig. [6l Figs. [Hand [7] for the 
dependencies of the diffusivities are the principal results of the current study. 

Of course, our planar triangle-like stars still represent a quite primitive system to 
mimic the non-ideal shape of real crowders in experimentally relevant setups. Future 
investigations including a three dimensional pyramid-like shape of crowders with longer 
polymeric arms will further elucidate the physical consequences of non-spherical and 
squishy crowders, and potentially exhibit additional unexpected behaviour. Moreover, 
not only the self-diffusion is to be studied but also the diffusion of tracer particles 
of various sizes and shapes in such crowded suspensions [56] as well as poly-disperse 
mixtures of crowders should be investigated. Some asymmetry may also be incorporated 
in the crowder shape. Recently, single particle tracking measurements allowed one to 
rationalise the translational and rotational diffusivities of micron-size symmetric and 
asymmetric boomerang-shaped particles in two dimensions |S7|. It was observed that 
the regimes of Brownian diffusion exist at short and long times while a coupling of D 
and Dr gave rise to subdiffusion at intermediate times. 
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